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Abstract 

We develop the formalism for canonical reduction of (1 + l)-dimensional gravity 
coupled with a set of point particles by eliminating constraints and imposing coordinate 
conditions. The formalism itself is quite analogous to the (3 + l)-dimensional case; 
however in (1 + 1) dimensions an auxiliary scalar field is shown to have an important 
role. The reduced Hamiltonian is expressed as a form of spatial integral of the second 
derivative of the scalar field. Since in (1 + 1) dimensions there exists no dynamical 
degree of freedom of the gravitational field (i.e. the transverse-traceless part of the 
metric tensor is zero), the reduced Hamiltonian is completely determined in terms of 
the particles' canonical variables (coordinates and momenta). The explicit form of the 
Hamiltonian is calculated both in post-linear and post-Newtonian approximations. 



1 Introduction 

Early versions of (1 + l)-dimensional gravity |l|, § have in recent years led to intensive 
study of a wide variety of such theories, in large part because problems in quantum gravity 
become much more mathematically tractable in this context f^. In contrast to (3 + 1)- 
dimensional gravity, the action integrals for these theories must incorporate some dynamics 
in the form of an auxiliary (or dilaton) field since the Einstein tensor is topologically trivial 
in two dimensions. Although most of these theories can be written in a generic form 
whose general properties can then be studied f^, 0] , there is often much to be learned by 
focussing on specific models within this general class 

One such theory has been extensively studied in many respects, including gravitational 



collapse, black holes, cosmological solutions and quantization |T0|, |TT], |T2|. Referred to 
as R = T theory, the specific form of the coupling of the auxiliary field \I/ to gravity is 
chosen so that it decouples from the classical field equations in such a way as to ensure that 
the evolution of the gravitational field is determined only by the matter stress-energy (and 



reciprocally) ^ |TOl- In this manner R = T theory captures in two spacetime dimensions 



the essence of classical general relativity (as opposed to classical scalar-tensor theories), and 
has (1 + 1) -dimensional analogs of many of its properties |T^, |TT|. Indeed, the theory can 



be understood as the D ^ 2 limit of general relativity |]13 . 

One of the fundamental problems of (l+l)-dimensional gravity is its relationship to 
Newtonian gravity. This is in general a problematic issue 0. For a system of particles the 
Hamiltonian in Newtonian gravity in two dimensions is 
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H = J27p^ + ^Gj2J2^am, \za-z,\ (1) 

a a h 

where ma, Za and Pa are the rest mass, the coordinate and the momentum of a-th particle, 
respectively, and G is the gravitational constant. We note that the R = T theory has been 
shown to have a Newtonian limit |ll|. However the dynamical role of the auxiliary field 



has not yet been fully analysed, a task which we investigate in this paper. 

More generally we shall, in the context of R = T theory, formulate a general framework 
for deriving a Hamiltonian for a system of particles, which coincides in a slow motion, weak 
field limit with the Hamiltonian (1). In addition to the properties stated above, an advantage 
of the R = T model is that it needs only one auxiliary field. As with the ADM formalism 
in (3 + l)-dimensional theory, we develop a canonical reduction by eliminating constraints 
and imposing coordinate conditions, which are quite analogous to the ADM conditions. In 
its final form the reduced Hamiltonian of the system is given as a spatial integral of the 
second derivative of the auxiliary scalar field Since in two dimensions there exists no 
transverse-traceless part of the metric tensor (which are the real dynamical variables of 
gravitational field), the scalar field \E' is given as a function of the dynamical variables [za, 
Pa) of the particles by solving the constraint equations. Then the Hamiltonian is completely 
determined in terms of the coordinates and momenta of the particles. 

The outline of our paper is as follows. Details of the reduction process are described 
in section 2, where much attention is devoted to the transformation of the total generator 
of the whole system and the choice of the coordinate conditions, from which the reduced 
Hamiltonian is defined. The consistency of the canonical reduction is proved in section 3. 
It guarantees that the canonical equations of motion given by the reduced Hamiltonian are 
identical with the original geodesic equations. The explicit form of the Hamiltonian is calcu- 
lated in section 4 in a post-linear approximation (a series expansion in the coupling constant 
G) and in section 5 in a post-Newtonian approximation. While the former approximation 
is appropriate for the analysis of fast motion, the latter is adequate for treating the slow 
motion, weak field case. Section 6 is reserved for concluding remarks. The role of the \E' field 
is discussed in the simplest case of a single static source. Finally, an appendix illustrates 
the relation between the equations of motion in our coordinate conditions and those in other 
conditions. 

2 Canonical formalism for particle dynamics 

The action integral for the gravitational field coupled with point particles is [|1^ 
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(2) 



where \l/ is the auxihary scalar field. Here Qfj^,^, g, R and are the metric tensor of spacetime, 
det^Qfj,^), the Ricci scalar and the proper time of a-th particle, respectively, and k = SvrG/c^. 
The symbol denotes the covariant derivative associated with g^iy. 
The field equations derived from the variations (5\E' and Sg^^i, are 

R - g^'V.V,^ = R- -^d^iy/^g^'d,^) = (3) 
^ V^^V,^ - ^g^^V^^Vx^ + (7^.V^Va^ - V^V,^ = kT^, (4) 



where 



T 



2 SC 



M 



y^nia dTa^=gf,agup-r--r^S^{x - Za{Ta)) , (5) 
„ J J-Q dTa dTa 



Cm being the matter Lagrangian given by the second term in the brackets on the right hand 
side of (2). The geodesic equation derived from the variation 6z^ is 



\9t.u{Za)-7^\ - -9uxA^a) -7^-7^ = ■ (6) 



The trace of Eq.(4) is 
which yields 



WV^^ = kT^^ (7) 

R = kT^^ . (8) 

Hence particle dynamics m R = T theory may be described in terms of the equations 
(6) and (8), both of which are independent of the scalar field. Particulate matter (in terms 
of the trace of its stress-energy tensor) generates spacetime curvature via (||); the effects of 
curvature act back upon matter via (|^). 

Note that all three components of the metric tensor cannot be determined from (^), 
since it is only one equation. The two extra degrees of freedom are related to the choice 
of coordinates. If the coordinate conditions are chosen to be independent of \E', equation 
(H) determines the metric tensor completely. More generally however, we need to know the 
scalar field ^ , through which the metric tensor is determined - it is this field that guarantees 
conservation of the stress-energy tensor via (^). So far little attention has been paid to the 
role of the scalar field ^ . 

To derive a Hamiltonian for a system of particles, we shall utilize the canonical formalism 
1^]. Writing 7 = gn,No = {—g^^)^^/'^,Ni = gio, the action (2) transforms to 



/ = Jdx^ i^Y.PaiaS{x - + 7r7 + H^ + NoR'' + NiR'^ (9) 
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where vr and 11 are conjugate momenta to 7 and \1/ respectively, and 

1 / ^' \' 



7' 1 



4/4^7 K a V 7 

(10) 



7 7 a 7 



Here and in the following we denote do by a symbol ( ' ) and di by a symbol ( ' ). 
Taking variations ^7, 611, 6Nq, 5Ni, 6"$, 611, 6za and 6pa of the action (0), we have 



(11) 



TT 



a 27^^ + m^ 



1 



+ iVi<'-^nvl/' + - + 5:^5(x-z,(x°)) 
I 7^ 7 a 7^ 

+ K „ V + N[- = 



2/^^77 



7 



7' 



7 - A/'o(2/t^77r - 2k^U) + Ni— - 2N[ = 



7 



1 



i?^ = 

1 



n + 5i(--iVin + — -iVo^' 

7 



«^V7 



iV^) = 



^ + No{2k^tt) - iVi(-^') = 



Pa + 



dNo Pi 
dza V 7 



— + mi- 



7 



iVo pI dj 



2jf + ml^'^'^ 



9Ni Pa Pa <97 ^ 

H ^Vi — 7- — — U 

dZa 7 7 "^^ra 

+ — = . 



2^ A^i 



4 - A^o 
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(12) 

(13) 

(14) 
(15) 

(16) 
(17) 



(18) 
(19) 



In the equations (18) and (19), all metric components (A^O) ^i, 7) are evaluated at the point 
X = Za and 

df _ df{x) 



dzn dx 



X — Za 



It can be shown that this set of equations is equivalent to the equations (3), (4) and (6). 
Here we show explicitly the equivalence of (18) and (19) to the geodesic equation (6). We 
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express the (/i = 1) component of (6) in terms of Nq, Ni and 7 and their derivatives. Using 
the relations 

1 

1 



we get 



d 



dt 




iVr 



IZa) 



+ - 



dZa 

'Nl 



. 7 7 dza 



dzr, 







(20) 



Ehminating pa from (18) and (19) yields an equation identical with (20). 

From the expression (9), the total generator (obtained from variations at the end point) 
is given by 



G = J dx < "^Pa^ix — Za)6Za + Tcdh — ^511 



(21) 



where the above form has been obtained by adding a total time derivative — 9o(n^) to the 
original action (9), and where the constraint equations (14) and (15) have been taken into 
account. In the above /i = 1 + 7; the advantage of using h will become clear later |T^ . 

We must now identify the dynamic and gauge character of the variables which appear 
in the generator (PT|). By considering the constraint equations (14) and (15), we see that 
the only linear terms there are and n' respectively. The equations may therefore be 

solved for these quantities in terms of the dynamical and gauge {i.e. co-ordinate) degrees of 
freedom. By writing 11 = ;^n", where 1/A is the inverse of the operator A = d'^/dx^ with 
appropriate boundary condition, we find that the generator becomes 



G 



27r' 



h-\ 




(22) 



where in obtaining this form, we have discarded surface terms. 

The form ( p2D of the generator is analogous what is obtained in (3 + 1) dimensions 
under an orthogonal decomposition of the hypersurface metric and its conjugate momentum 



1^, . We therefore propose adopting the coordinate conditions 



X 



t 



2A 



-h' 




(23) 
(24) 



which will then allow us to identify the Hamiltonian and momentum densities as the re- 
spective coefficients of 6t and 6x in the canonical form of the generator. Taking a spatial 



5 



derivative of these coordinate conditions yields their differential form, which does not ex- 
plicitly depend upon the coordinates {t,x). This leads to the choices 



A 



h = 2 

n' = 



7 = 1 

n = 



(25) 
(26) 



which may then be inserted in equations (12) and (19) to solve for the relevant physical 
degrees of freedom. 

Specifically, we may insert these choices into (14) and (15) to solve for the Hamiltonian 
and momentum densities. We then find that the generator (^) may be expressed in the 
canonical form 



G = J dx i^Pa6{x - Za)6Za - %f,6x^^ 



where 



T 

-^00 



n = 

271' 



(27) 

(28) 
(29) 



Thus, Ti is the Hamiltonian density of the system and 7^^ is the momentum density. Note 
that since the constraints (14) and (15) have already been imposed, %^ is expressed in terms 
of the canonical variables Za and pa as '7o^(x, Za,Pa), by solving these constraints. With the 
coordinate choices (25) and (26), they lead to 



2ti' + ^pj{x - Za) 







(30) 



(31) 



Following a similar procedure as that used to transform the generator G, we rewrite the 
action integral (9) as 



j (fx l^PaZa5{ 



27r' 



Iff 



h-l A 



h-l 



(32) 



where we have discarded surface terms. Eliminating the constraints (14) and (15) and 
imposing the coordinate conditions (23) and (24), the action integral reduces to 



Ir = J dx^ \^PaZa5{ 



X - Za ) - Hj . 

Thus the reduced Hamiltonian for the system of particles is 

1 



H 



dxH 



K 



dxMf 



(33) 



(34) 
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where \E' is a function of Za and Pa and is determined by solving the constraint equations (30) 
and (31). This expression is analogous to the reduced Hamiltonian in (3 + 1) dimensional 
general relativity [T^, |l^, |l^ . 



We pause to comment on the relationship between the integral forms (23) and (24) of the 
coordinate conditions and the differential forms (25) and (26). These two forms are equivalent 
only when one retains the appropriate boundary conditions for the integral operator 1/A. 
A proper treatment [p!7| , |18| entails the insertion of a regulator exp(— a|a;|) in the left-hand 



sides of (23) and (24), yielding 

^:=2e(-"l^l)(l-a|x|)- 

and 




2e(-°l^l) -2|x|ae(-"l^l) -ly/^Kj 

in place of (25) and (26) respectively. We must then insert these quantities in (12) - (19), 
taking the limit a — > at the end of the calculation (implicitly assuming a\x\ < 1). This 
turns out to be equivalent to inserting (25) and (26) in these equations. The action (33) is 
recovered by a similar limiting procedure. 

This situation is analogous with that in (3+1) dimensions dimensions as discussed by 
ADM in ig (cf. Eqs. (80), (81), (84) and (85) in section 6). We shall discuss these points 
further in section 6. Note that the coordinate conditions in integral form are needed only at 
the stage of transforming the generator and defining the Hamiltonian density. After fixing 
the formalism, we need only the differential forms to solve the constraints and to evaluate 
the Hamiltonian, a problem which is treated in the following sections. 



3 Consistency of the canonical reduction 



Following a procedure similar to the original ADM argument [^, IT^j , we shall demonstrate 
here that the canonical equations of motion derived from the reduced Hamiltonian (34) are 
identical with equations (18) and (19), which are in turn equivalent to the original geodesic 
equation (6). 

We start from the action integral (9). The variations with respect to Pa and Za lead to 
the equations of motion 

'''yNJ,)'-^ (35) 

These equations are identical with equations (18) and (19). Our purpose is to prove that 
these equations lead, when they are combined with the constraints (14) and (15), to the 
equations in the reduced formalism. 
Defining 

Po{x) = --A^(x) (37) 

K> 

Pi{x) = 27r'(x) (38) 
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and imposing the coordinate conditions, the action (p2D may be expressed as 



I = j ^^PaZa^i 



(39) 



where we have not yet imposed the constraints = 0. By solving the constraint equations, 
we get 

Poix) = %oix,Za,Pa)=n (40) 
Pl(x) = Toi{x,Za,Pa) (41) 

where Tq^{x, Za,Pa) was described in the previous section. 

Expanding in a functional Taylor series about the point = Tq^ gives 

'dR^ix) 



R'^ix) = J d^y [P,{y)-%,{y,za,pa)] 



SPuiy) 



+ 



(42) 



Jp=r 



where the coordinate conditions have already been imposed. Substituting this expansion 
into the right hand side of (39) and taking the variation of the action / with respect to P'^, 
yields 

'SR^^iy) 



-6,0 + J d^yN^iy) 



6PJx) 



+ ■ 







(43) 



J p=r 



where the terms represented by ■ ■ ■ contain [P,^ — Tq^] as a factor. 

Now we need the relation which is valid after imposing the constraints. Then, by requiring 
R'^ = 0, we have 

'SR^y) 



6,0 + J d^yN^iy) 
Insertion of (42) into (35) leads to 



6P,{x) 



. 



(44) 



J p=r 



Za 



d'yN^iy) / d^x 



dTo 



dpait) 



6R^{y) 
6PJx) 



p=r 



+ {Pu{x) -To,{x,Za,Pa)) 



d 



dpa{x) 



6Rny) 
6P,{x) 



+ 



J p=r 



Imposing the constraint equations i?^ = yields for (45) 



d X 



dTo 



dpa{t) 



d'yN,{y) 



6R^{y) 
6PJx) 



(45) 



(46) 



J p=r 



As a consequence of the relation (44), (46) becomes 



d X 



dTo 



00 



dpa 



OH 

dpa 

In the same way we can show that the equations of motion (36) is identical with 

dH 



Pa 



dZa 



(47) 



(48) 



Thus the consistency of the reduced canonical formalism is proved. 
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4 Post-linear approximation 



For a direct calculation of the reduced Hamiltonian (34), we have to solve the constraint 
equations (30) and (31). However, it seems quite a difficult task to get an exact solution 
except in the case of a single static source. So we need some approximation method. 

In this section we apply an iterative method successively to the integrand of the expression 
(34) with the use of the equations (30) and (31). This provides us with an expansion of (34) 
in powers of k, which we refer to as the post-linear approximation. 

First, substituting A^' given by the equation (30) into the integrand and performing a 
partial integration, we have 



where x is defined by x' = vr and 



(49) 



1 1 °° 

-— + KXX 
4k J _o< 



(50) 



Noting that ^ is of order k, we substitute (30) for A\E' and (31) for Ax into the integrand 
of the second term on the right hand side of (49). Iterating twice in terms of k yields the 
expression 

+ ^ EP'^X(^a) |l - ^*(^a)} + y E + X{Za)x{Za) + ^2 



where 



So 



K 



X 



32k 



(52) 



The last term in ( |5lD contributes at order k^ . This iteration can be continued successively 
and corresponds to a series expansion in k. 

In contrast to the (3 + l)-dimensional situation, we encounter a subtle problem when we 
try to extract an explicit form of the Hamiltonian under this iteration scheme. In (1 + 1) 
dimensions the surface terms such as 5*1 and 5*2 arising in the process of the calculation do not 
necessarily vanish and make the Hamiltonian indefinite. This is because the dimensionless 
potential Gmr/c^ becomes infinite at spatial infinity, in contrast to (3 + 1) dimensions where 
the corresponding quantity Gm/rc^ vanishes at spatial infinity thereby assuring that the 
associated surface terms vanish. It is therefore an important task to find a solution of \E' and 
X which makes the surface term vanish at a given order in k. This is related to the choice 
of boundary condition. 



9 



From the expressions of and S2, we can infer a simple choice of the boundary condition 
as follows. 

For f{x) = *2 - 4:K,^x^, 

f{x) = and f'{x) — in a region | x \»\ Za \ for all a. (53) 

It is easily checked that under this boundary condition, the surface terms and 5*2 exactly 
vanish, since they are proportional to /' and (4 — ^)/' + /^f' respectively. 

Now let us try to obtain the Hamiltonian up to (the 2nd-post- linear approximation) . 
First we expand and x in a power series in k 



X^'^ + KX^^) + ■ • • . 
Substituting this expansion into the equations (30) and (31), we get 



X 



A^W = -J2^pl + mlS{x-Za) 

a 

^ a 

Ax^'^ = . 

The solutions which satisfy the boundary condition (53) are 

= -\ J2 [para + e^pl + ml {x - Za)| 



(54) 
(55) 



(56) 

(57) 
(58) 
(59) 

(60) 
(61) 



^(2) 



1/1^2 

"2 u 

2 



a 



J2 VPI + mlra + € Y.Pa{x - 



\IpI + ^ab + e YPb{Za - Zb) 
b b 



-J2Po-^'^ 



J2Pb1"ab + e J2yPb + 1TT'b i^a - Zb) 
. b b 



e (\ 



2 V4 



NPI ■^'ml{x- Za) 



Y^Pl + ^rab + e YP^iza-Zb) 

Lb b 



J2Pa{x - Za) 



YPbfab + e YyPb + ^b W - Zb) 
Lb 6 



where ra=\x — Za\- 
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In these solutions we introduced a constant of integration e, satifying = 1. We have 
two types of solutions corresponding to e = 1 and e = — 1, which are related to each other 
under time reversal. This e factor guarantees the invariance of the whole theory under the 
time reversal. 

The boundary condition is checked as follows. Up to order we have 

,2 



T 



EE 

a b 



Pi + ^ly Pb + - PaPb 



+j\/pI + KY1 



Pb + ^bc - PbBbc 



ran - {X- Za){x - Zb)] 



K 



+^ YEE|vPa + "^a (^1- gE 



p1 + rnlrc + e Pc{x — Zc) 

-^Pa E \/Pb + '^b ^bc - PbBbc \ [raix - Zb) - {x - Za)rb] 



where 



Abe = ^Jpl + ml Tbc + e Pc{zb - Zc) 
Bbc = Pcnc + e VP? + ml{zb - Zc) 



Since both [raTb — {x — Za){x — Zb)] and [ra{x — Zb) — {x — Za)rb] vanish in a region | x \»\ 
Za \,\ Zb \, the boundary condition is satisfied. 

The 2nd-post-linear Hamiltonian is therefore unambiguously determined to be 

H = E\/pl + <{^- l^i^a) + ^^{Za)^{Za)] 

+ EP«X(^a) {l - l^(^a)} + y E \JpI + ^1 X{Za)x{Za) 
= E V^a + + ^ E E (JpI + K^JPb + ^6 - PaPb) Tab 
+ Y E E (\/Pa + TT^l Pb - Pa\/Pb+ ^b) {Za - Zb) 



a b 



4 V4 

EPa 



EV^ + \ j2Pbfab + e E \/Pb+ ^bi^a - Zb) 



a Lb 



a b 



J2Pbfab + e E 'JPb+'^bi^a - Zb) E + ^ac + ^ EPc(^a " ^c) 

. b b lie 

+ EE \/pI + ^l\IPb + rnl Tab - e Pa\IPb + ml{za - Zb) 

X E + nc + e ^Pc{Zb - Zc) 



a b 



-EE \IpI + Pb^o-b - ^ PaPb{Za - Zb) Y^PcTbc + e X] V^Pc + "^c(^6 " 
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(64) 



As shown in the above development, the K-expansion is the successive approximation 

in the background of Minkowskian space-time and the terms in each order of k preserve 
relativistic forms. This approximation is therfore appropriate for describing relativistic fast- 
motion of the particles and can be carried out to any desired order of n. 



5 Post-Newtonian Hamiltonian and the redefinition of 
canonical variables 

To compare R — T theory with Newtonian gravity, we need an approximation method 
applicable to a slow motion in a weak field. It is provided by the so-called c"^ expansion. 

All terms appropriate to the post-Newtonian approximation (up to the order of c~^) are 
included in the post-linear Hamiltonian (64). Noting that both p^/m^ and ■\/k are of the 
order of c~^, we find from (64) 

H = Y.rna + J2w^+oJ2J2^ambrab + —J2J2i^aPb-mbPa){za-Zb) 

a a ^ a b ^ a b 



4 ^ 2 ^ 



+ \ (^) H n^ambmcTabrac + ^ (^jj mm mambniciZa - Zb){Za - Zc) . 

(65) 

The corresponding solutions of and x 

^ ^'^^ Y,ma{x-Za)] - fj) J2J2'^o.mbrarab (66) 



2\AJ W V4 



a b 



^ a ^ a 

The canonical equations of motion yield 



dH 

Za r-i 

OPa 




_ 

rria 


— }_^mb{za- 

^ b 


dH 

OZa 






drab ' 
dza 



3 

Pa I ^ Sr^ Pa ^ sr^ /no\ 



ab 
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~^ (^) Z!Zl"^«"^6"^c(^a - ^b) ■ (69) 



By eliminating pa in the equations (68) and (69), we get the equations of motion to second 
order 

rriaZa 

drab f./. n2 n2~1 drab 



4 V dza 8^ ^ ^ dza 



kV^^ fdVar. drbc\ 1 //^A^V^V^ f ^Tab , ^r^fc 



n \ \ ^ \ ^ I ' ac be \ \ \ ^ \ ^ / ab 



1 ^ 

+ 2(4) J2J2'^o.'^bmc{za- Zb) (70) 
6 c 



We pause to show that these canonical equations of motion arc equivalent to the geodesic 
equation (6). To evaluate the geodesic equation, we have to know all components of the 
metric tensor, of which gu = 7 is fixed and Nq and Ni components are determined by 
combining the time derivatives of the coordinate conditions (25) and (26) with the equations 
(13) and (16). These equations are 

Kx'No + N[ = (71) 
^iQiVo^' + TV^) =0 . (72) 

The solutions in the post-Newtonian approximation are 

K y-^ eK ^ K y-^ Pa 



+ |(^) (j^^aix-Za)^ ^^(l) JLJL^o-^bTarab (73) 

^1 = T 5Z ^«(^ - + 7 E Par a ■ (74) 

The boundary conditions in the previous section are incorporated in these solutions 
through (12) and (17), which serve as consistency equations. 

Since we have proved in section 2 the equivalence of the geodesic equation (6) to a set of 
equations (18) and (19), we have only to evaluate the latter equations. Under the coordinate 
conditions (25) and (26), the equations (18) and (19) become 

= NoiZa) - N,{Za) (75) 



Pa - yPa + rriag^^ + Pa (76) 
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which lead to the equations (68) and (69), after the insertion of the solutions (73) and (74) 
and expanding in powers of c~^. 

The Hamiltonian (65) contains a term proportional to c^^, namely, 



a b 



— J2 Y^if^aPb - mbPa)iZa - Zb) . 



The appearance of such a term may seem unnatural, because in (3 + 1) dimensions terms 
in odd powers of , for example terms, are related with gravitational radiation; yet 
there exists no graviton degree of freedom in (l+l)-dimensions. However, by a suitable 
redefinition of canonical variables 

Za ^ Za Za 

(77) 



Pa > Pa = Pa - ^ambiZa - Zb) 

^ b 



{71 



this term can be eliminated. Under this redefinition the "Poisson brackets" among the 
canonical variables are kept unchanged. The Hamiltonian expressed in terms of the redefined 
canonical variables is then 



H = i:ma + T.^ + li:i:mambrab-T.A + ^j:T.ma^rab 

a a O a b a ^"^a ^ a b 

~^ mi^a^^'^afe + 7 fy) J2YY'^'^'^b'mc{fabrac- {Za- Zb){Za- Zc)} (79) 

° a b ^ ^^-^ a b c 

It is straightforward to show that the equations of motion derived from the Hamiltonian (79) 
are identical with those derived from the original variables. 



6 Discussion 

We have performed the canonical reduction of a set of point particles coupled to a gravita- 
tional field in (1 + 1) dimensions in the context of the R = T theory, and have shown how to 
carry out series expansions of such equations in powers of k (the post-linear approximation) 
and powers of (the post-Newtonian approximation). We have obtained explicit expres- 
sions for the Hamiltonian to order and respectively, along with an expansion of the 
geodesic equation to order c""^. 

We have considered in this reduction the dynamical role of the auxiliary scalar field 
Despite the decoupling of \1/ from the classical gravity /matter field equations, the Hamil- 
tonian of the system is defined in terms of One of the essential points in this analysis 
is the choice of the coordinate conditions, which we chose in the forms of (23) and (24) 
through the transformation of the generator. We note that our coordinate conditions are 
quite analogous to the ADM coordinate conditions in (3+l)-dimensional gravity WM. These 
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coordinate conditions (frequently used for particle dynamics [|T^, |l^) are 




80) 



81) 



where 




9ii A 9ij,ij 



1 

A 




82) 



83) 



and 71^^ is a conjugate field to Qij, A being the Laplacian operator in three dimensions. 
The differential forms of the conditions (80) and (81) are 



^9ij,j ^9jk,ijk ^^9jj,i — 

TX = U . 



84) 
85) 



The first condition (84) relates space-derivatives of spatial components gij of the metric 
tensor. In (1 + 1) dimensions there is only one spatial component 7 = gu] from (23) the 
condition analogous to (84) is 9i7 = 0. Therefore we chose the condition (25). Since vr*-' 
in the second condition (85) is the conjugate momenta to gij, the corresponding condition 
in (l-l-l)-dimension might at first appear to be vr = 0. However, the form of the generator 
(21) indicates that when the variation 5h is related to one of the coordinate conditions, 
the variation of the conjugate vr can no longer be exploited for the other condition. The 
candidates would be 511 or 5^ . Judging from the fact that vr*-' is essentially a time-derivative 
of gij and 11 is also a time-derivative of 7, we have no choice to take other than 5TI. 

We turn now to a discussion of the consistency of the full set of equations (12) - (19). 
As discussed in previous sections, to determine the Hamiltonian for a system of particles, 
we have only to solve the constraint equations (14) and (15) for \1/ and vr. All informations 
on the dynamics of particles are included in the constraint equations. Equations (13) and 
(16) are used to determine A^o and Ni, and equations (18) and (19) lead to the equations 
of motion. The remaining two equations (12) and (17) are not necessary in describing the 
dynamics of the particles. However they are necessary in providing a full description of the 
gravitational field. 

To check the consistency of the whole formalism, consider substituting the post- linear 
solutions of \E',7r, A'o and A^i into (12) and (17). After some calculation, it is straightforward 
to see that (12) is consistently satisfied but (17) is not. The reason is as follows. In the 
general solution to the constraint equation (14) (or (30)) \I' may contain in general an 
a;-independent function f{t). Since all other equations except (17) contain only spatial 
derivatives of f(t) does not contribute to either the Hamiltonian or to the equations 
of motion. The function f{t) is necessary only for the consistency of (17). The boundary 
condition and the consideration on the surface term in the section 3 can be applied with \l/ 
replaced by (\I' — /(t)), the latter function having no explicit dependence on t. After lengthy 
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and complicated calculation of the equation (17), we can determine the explicit form of f{t) 
to the order k^. 

To illustrate this situation, let us consider a single static source at the origin. In this 
case the fundamental equations are 

= KTT^ - ^ {<i/f + MS{x) (86) 

tt' = (87) 

K7rNo + N[ = (88) 

ai QiVo*' + iVo) = (89) 



7r + A^o 



3fi; 9 1 ,^,x2 



+ Nitt' + T^A^o*' + N[Tr^O (90) 



4f + 2fi:iVo7r - iVi^' = . (91) 



Solutions to the constraint equations (86) and (87), which satisfy the same boundary condi- 
tion in the section 3, are 

* = -^l^l+/(^) (92) 
^ = X' = - (x = -'-fx) . (93) 

Here we included f{t) in ^. 

The solutions of Nq and A^i are 

No = e^l-l (94) 
N, = ,^(e^N_i) . (95) 

Equation (90) is satisfied, whereas (91) leads to 



/(^) = ^ ■ (96) 



Then the solution of ^ becomes 



* = I X I H — — i . (97) 



This shows that even for a static source the auxiliary field is not static. 
The Hamiltonian is 

= dxA^ 

K 



M . (98) 
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Since 



^00 = -iVo' + ^i 
< -1 

this solution has no event horizon. 

The canonical formalism we have derived can be utilized in a number of ways. The 
most obvious of these is to obtain an explicit solution for the A^-body problem in (1 + 1) 
dimensions (the preceding discussion illustrates the solution when = 1). In the large- 
limit with fixed proper distance between the most widely separated particles one might 
expect to recover the fluid collapse problem studied in ref. [0]. More generally, one could 



study the iV-body problem where gravity is coupled to additional matter fields. These all 
remain interesting subjects for future investigation. 



Appendix 

Relation between the equations of motion in the 
coordinate conditions 7 = 1, 11 = and those in the 

conditions ^00 = -^fiS^oi = 

One of the most popular coordinate conditions in (1 + 1) dimensions is the one for which 
the metric tensor has the form 

= ( T a° . ) • (89) 

The difference between this condition and the condition we adopted in this paper is that 
the former is independent of the auxiliary scalar field \E' but the latter is not. This \E'- 
dependence of the coordinate condition played an important role in the canonical reduction 
and the general expression of the Hamiltonian. 

In this appendix we shall compare the equations of motion in these two coordinate con- 
ditions. Expanding a and the energy-momentum tensor (5) in a power series of we 
obtain 

« = l + 20-(7SJ + --- (100) 
Too = E / o?l^?>{x-Za) 

= Y.^a[l+^<P+\{Zaf + --^S{x-Za) (101) 

X ^ rriaOL I Za r-, s 

Toi = . b[X-Za) 

a y^l-a~2(ia)2 

= Y.^a{-Za + ---}5ix-Za) (102) 
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Til = l" " S{x-Za) 



= Y,ma{z,y{l + ---}6{x-z,) (103) 

a 

where the number in the upper parenthesis denotes the order of and = — ^^'oo''- 
The post-Newtonian equation of motion of the reference |1T| reads 

mi'zi = —mi(f)' — mi^p' — 4mi(p(p' + Srriicpzi + 3mi(f)'{ziY . (104) 

Here (p and ip satisfy the equations 

r — 2 00 

= -Xl"^a5(a; - ^a) (105) 

^ a 

r = -0-2(0')' + |(r°°('^-T"(2^) 

= -4>-2{<pr-'^Y.^al^^Y.^trat + l{zaY'^5{x-Za) . (106) 



The solutions are 



AC 

a 

^ a ^ a 

- ) rnanibran + 2 ( - ) '^arribrahra 

a b a b 

-LO a 6 O a 

The equation of motion (104) for a system of two particles becomes 

miZi = --mim2-z —mim2{mi + m2){zi - Z2) 

4 ozi o 

+ ^mim2{6(ii)2-6iii2 + 3(i2)'}^ • (109) 

To investigate the relation between this equation of motion and the equation obtained in 
our canonical formalism, we consider coordinate transformations x'^ — > xj^, which connect 
two forms of the line element : 

ds^ = -adt^ + -dx^ (110) 

a 



{N^ - Nf)dr + 2Nidtdx + dxf , (111) 
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The explicit form of the transformations is given by 

i = t-'-^J2^a{x-Z,f + --- (112) 

O a 

X = X+-^ma{x-Za)\x-Za\^ ■ (113) 

O a 

For a system of two particles the above transformation (113) leads to 

Zl = ^1 + ^rn2(ii - I2) Ui - ^2 I H 

^2 = ^2 - ^"^i(li - I2) Ui - ^2 I H 

^1-2:2 = (^1 - I2) |l + ^("^1 + "^2) Ui -I2 I H 1 

^12 = £12 |l + ^("^1 + "^2)1:12 H 1 • 

Since the second term of the right hand side of the transformation (112) contributes to f as 
a correction of the order c~^, in the post-Newtonian approximation we can set t — t. Then 
we have 

Zl = ii + |rn2 |(ii -^2) Ui -^2 1 +(^1 -^2)(ii -^2)^^! + • • • 

= ii + ^"^2(ii - 'Z2) Ui - 12 1 H — 
Zl = ii + 7^2 1 (li - 'Z2) Ui - ^2 1 +(ii - ^2)^^ 1 H — 



^2:12 



5ri2 d 

|1 - -(mi + m2)ri2 + " " 7 |^ + 4 (^1 + ^2)^112 + 
+ 0(k^) . 



1 + - {rni+ 7712)112 + 
o 



9^.1 

By the use of these relations the equation of motion (109) is transformed to 

771ili = -T"^l"^2^=^ + - {^{Zif - 2iii2 + (^2)^ ^^-T^"^1^2(?77l + "^2)(^l-^2) (114) 

4 azi 8 ' ' J c/^i 16 

This equation is identical with the equation of motion (70) we obtained in the canonical 
formalism. 
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